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The redshift dependence of the observed cluster abundance due to the Sunyaev-Zeldovich 
effect (SZE) is studied in various cosmological models, including flat and open homogeneous 
(CDM) models and an inhomogeneous model with a large-scale local void. The Press- 
Schechter formalism is used to derive the abundance at epochs in the range < z < 2, 

■ and the cluster mass limit Mij m is obtained from a flux limit for SZE. It is shown that 
£SJ ' SZE is useful for constraining the cosmological model parameters, and the abundance in the 

inhomogeneous model may be comparable with that in the low-density homogeneous models. 
The significance of relative difference of abundances in the various models is discussed. 

> ' 
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. §1. Introduction 

■ The abundance of clusters has been extensively studied by many people in X-ray 
surveys to constrain the cosmological parameters: e.g., Bahcall and Fan, 1 ) Viana 
and Liddle, 2 ) Eke, Cole and Frenk, 3 ) and Kitayama and Suto. 4 ^ 6 -* The cluster 

^ ■ abundance found in the submm survey based on the Sunyaev-Zeldovich effect (SZE) 

has also been studied by Haiman, Mohr and Holder, 7 ) Holder, Haiman and Mohr, 8 )' 9 ) 
Kitayama, Sasaki and Suto, 10 ) and Fan and Chiueh. 11 ) This SZE survey is expected 
» ■ to be most important in clarifying the evolution of clusters and constraining the 

\ cosmological parameters. 

In this paper we study the cluster abundance obtained in the SZE survey, whose 
observational conditions were set to correspond to the interferometric arrays in the 
AMIBA project. The cluster mass limit M\[ m for deriving the abundance is deter- 
mined using the expression for the flux S v given by Kitayama and Suto, 5 ) which 
is different from that used by Fan and Chiueh. 11 ) This leads to a difference in the 
behavior of the resulting cluster abundances from that reported in their paper. 

In §2, we describe the formulation for deriving the cluster abundance in the 
SZE survey. In §3 we present the results in various cosmological models. Here, we 
first consider four representative homogeneous cosmological models: LCDM with 
(/3b, A ) = (0.3, 0.7), h = 0.7 and r = 0.25, OCDM with (0.3, 0), h = 0.7 and r = 
0.25, SCDM with (1.0,0), h = 0.5 and r = 0.5, and rCDM with with (1.0,0), h = 
0.5 and r = 0.25, where r is the CDM shape parameter and Hq = lOO/i km s~ x 
Mpc -1 . 

At present, most cosmological observations, including SDSS, 12 ) high-redshift su- 
pernovae 13 ) -16 ) and WMAP, 17 )' 18 ) support a flat homogeneous model with nonzero 
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cosmological constant. However, the presently observed values of the Hubble con- 
stant seem to be non- uniform, 19 )~ 24 ) though they include rather large uncertainties; 
that is, the local median value seems to be larger than that of the global median 
value by a factor approximately equal to 1.2. If the non- uniformity of the Hubble 
constant is found to be real, we may have to use inhomogeneous models to describe 
the cosmological observations. In the magnitude-redshift diagram of supernovae of 
type la, moreover, the point representing the recent data 26 ) with z = 1.7 deviates 
from the curve predicted by the concordant model with nonzero cosmological con- 
stant. If additional supernovae with z > 1.5 exhibiting a similar trend are obtained, 
models different from the concordant model may be needed in order to account for 
the data. Taking this situation into consideration, we consider here an inhomoge- 
neous model with a large-scale local void as a representative inhomogeneous model, 
in which there are inner and outer homogeneous regions (I and II) and a spherical 
boundary. 27 )~ 30 ) 

The objectives motivating our study of this inhomogeneous model are discussed 
in previous papers. 31 )' 32 ) The first of these objectives is the introduction of the inho- 
mogeneity of the Hubble constant, and the second is the explanation of the acceler- 
ating behavior of SNIa in the absense of a cosmological constant. The cosmological 
parameters in the two regions are (1?q,Aq) = (0.3,0), r 1 = 0.25, Hq = 70 km s _1 
Mpc- 1 in the inner region (I) and (f$, Aj, 1 ) = (1.0, 0), T n = 0.5, H$ = Hi x 0.82 in 
the outer region (II). It is assumed that our observer is at the center for simplicity, 
and the spherical boundary corresponds to a redshift Z\ = 0.067. For the above 
set of parameters, the observed magnitude-redshift relation of SNIa, including the 
above high redshift data, was reproduced in the inhomogeneous model. 30 ) The con- 
sistency of spherically symmetric inhomogeneous models with the supernova data 
was recently discussed and examined also by Iguchi et al. 33 ) In §3 we discuss the 
difference among the abundances in these various models. In §4 we have concluding 
remarks. 



§2. Number density and SZE 



2.1. The observed number of clusters 

The comoving number density of clusters of mass M with width dM is 

n(M)dM = (I) V2 ^^#^ exp (2-1) 

following the Press-Schechter formalism, where po is the present mass density of the 
universe, 5 c (z) is the linear density threshold for collapse at redshift z, and oo is the 
rms linear density perturbation on the scale corresponding to M. The expressions 
for S c (z) in homogeneous models are given in Kitayama and Suto's paper 5 ) (cf. their 
Appendix A). 

The differential number of clusters is expressed as 



dN dV 



dzdfi dzdfi J M 



[ n(M)dM, (2-2) 
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^o = ct 8 (z)( — j— ) , (2-4) 



where dfi is the solid angle element, dV is the comoving volume element, and M\{ m is 
the lower limit of the observed cluster mass, which is discussed in the next subsection. 
Here, the comoving volume element dV is given by 

dV = ct [d A (z)] 2 d[2(l + zf, (2-3) 

where cZa is the angular diameter distance. Following Viana and Liddle, 2 ) we assume 
(To in the form 

R \ -7(H) 

,8/i^Yipc/ 

where 7(12) = (0.3r + 0.2)[2.92 + log 10 (i2/8/i _1 Mpc)], with M = fvrpo^ 3 , and we 
have 

o- 8 (z) = a s (0)g((2(z), A(z))/[(1 + z) g(f2 , A )], (2-5) 

where g(f2, A) is an approximate factor representing the growth of linear density 
perturbations (Carroll, Press and Turner 34 )), given by 

g(f2, A) = ^n/[f2 4 / 7 - A + (1 + fl/2)(l + A/70)], (2-6) 

and f2(z) and X(z) are ft and A at the epoch z. 

For homogeneous models the observed values of 08(0) have been derived by 
several groups: for instance, cr 8 (0)^- 45 = 0.53 ±0.05 for A = and <r 8 (0)^g- 53 = 
0.53 ± 0.05 for f2 + A = 1 by Pen 35 ) and Eke et al., 3 ) cr 8 (0)^- 60 = 0.50 ± 0.04 
by Pierpaoli et al., 36 ) <7 8 (0)l2g- 6 = 0.345 ± 0.05 by Fisher et al., 37 ) <7 8 (0)J3§- 5 = 
0.33 ± 0.03 by Bahcall et al., 38 ) cj 8 (0)^- 48 ~ a27 ^ = 0.38 by Viana, Nichol and 
Liddle 39 ) and <r 8 (0)/2g- 5 = 0.48 ± 0.12 by Spergel et al. 18 ) for the WMAP data. 
Here we adopt the following two sets of values: 

o- 8 (0)flg = 0.4 and 0.5, (2-7) 

where pis 0.45 and 0.53 for models with Ao = and flat models with J?o + Ao = 1, 
respectively. 

2.2. The lower mass limit M\- im 

When photons pass through a cluster of hot electrons, a temperature decrease 
results and the black-body spectrum is distorted due to inverse Compton scattering 
as 

AT/Tqmb = g{x)y, (2-8) 
where the Compton y-parameter is 



V 



^ )s T^m- 4 - (240) 

Here, x = hu/(kTcMB),whereu is the CMB photon frequency, n e is the electron 
number density, <tt is the Thomson cross section, T gas is the temperature of the 
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cluster gas, and the integration is along the line of sight. If T gas 3> Tcmb> the flux 
of CMB photons changes from S™ B (= (2hv z j <?) j (e x - 1)) to 

Su = S™ B ^-g(x)Y, (2-11) 
e x — 1 

where 



Y = [d A ]- 2 J ydA, (2-12) 



with dA the element of the projected area of a cluster. For v = 219 GHz, we have 
g(x) = 0. For v < 219 GHz, g(x) < 0, and for v > 219 GHz, g(x) > In the Array 
for Microwave Background (AMIBA) project (Fan and Chiueh 11 )), the parameters 
are v = 90 GHz and x ~ 1.58, which we use in the following. For an isothermal 
cluster with constant gas mass fraction, we have 

Y = -^^[d A ]- 2 f 1CM (l+X)kT gas M, (2-13) 
2m e m p 

where m p is the photon mass, X is the hydrogen mass fraction, and /icm = ^b/^o 
for the present baryon density parameter Q-q. For an isothermal gas, T gas is related 
to the total cluster mass M by 

as shown by Kitamura, Sasaki and Suto, 10 ^ and the expression for p v i r {z) (which is 
the ratio of the mean density of the virialized cluster to the mean density of the 
universe at each epoch) is given in the paper of Kitamura and Suto 5 ) (cf. their 
Appendix A). Here, z is the redshift for cluster formation in principle, but it is 
regarded here as the redshift at the epoch in which the cluster exists. Then, the 
total flux of the cluster is 



S - 25 5Ml + z)a(*)i±^ x ^ Mz)^ ( Pvh{z) y /3 ( M ,5/3 
b u -2b.bh{l + z)g{x) 2 x ^ /3 [ c/Hq \ UsttO U^-imJ mJy ' 

(2-15) 

where g(x) is given by g(x) = x 4 e x (e x — 1) 2 g(x). This expression is different from 
that given in Fan and Chiueh's paper, based on Eke, Cole and Frenk, 3 ) mainly with 
respect to Q$. This difference seems to result from the difference in treating the gas 
mass fraction as given by Q-q/Qq or as fixed. 



§3. Cluster abundance in various models 



In this section we assume that the limiting flux (S v )n m is 6.2 mJy, corresponding 
to the AMIBA design, and calculate the limiting mass M\\ m using Eq. (|2-15|) . Then, 
by integrating Eq. ()2-2|) . the differential number density of SZE clusters, dN/ (dzdf2), 
is obtained. In the following, we carry on this calculation for various models men- 
tioned in §1, namely, the homogeneous models (LCDM, OCDM, SCDM and rCDM) 
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and an inhomogeneous model with two homogeneous regions I and II, which cor- 
respond locally to OCDM and SCDM, respectively. The four homogeneous models 
correspond to (i?o,Ao) = (0.3, 0.7), (0.3, 0), (1.0, 0) and (1.0,0), respectively, and 
r = 0.25, except for SCDM with r = 0.5, where r is the CDM shape parameter. 

3.1. Homogeneous models 

In Figs. 1 and 2 the types of behavior exhibited by dN/(dzdQ) in the four 
homogeneous models are shown for the two sets (<78(0)I?q = 0.5 and 0.4) given in 
Eq. (|2-7|) . Here, we use units of deg~ 2 and set Q-^h 2 = 1.70 x 10 -2 . It is found 
that (1) the peaks in the low-density models are at epochs z = 0.17 and 0.07 for 
C8(0)J2g = 0.5 and 0.4, respectively, and the number density decreases as a function 
of z more slowly for larger erg(O), and (2) for z = 0.2, the number density in LCDM 
is larger by factors > 10 than that in SCDM in both cases of o"s(0). This trend of 
latter behavior is the same as that seen in Fig. 3 of the paper of Holder et al., 8 ) but it 
is different from the behavior displayed in Fig. 1 of the paper of Fan and Chiueh. 11 ) 
This difference comes mainly from that of the factor Qq used in the expressions of 
S u in the two approaches. 

Next, we consider the ratio r of the number [iV(< 0.5)] of clusters with z < 0.5 
to the number [iV(> 1.0)] of clusters with z > 1.0 in various models, following 
Fan and Chiueh. In Table I, the ratios for OCDM and LCDM are listed, while the 
ratio for SCDM is omitted, because it is very large compared with that for the other 
two models. It is found that the ratio r is sensitive to the value of crs(0) as well 
as the model parameters. Accordingly, the observation of r may yield a significant 
constraint on them. 

Table I. The ratio r = N(< 0.5) /N(> 1.0) in various models for two values of <ts(0)(J?o) p , where p 
is 0.45 and 0.53 for models with Ao = and flat models with Qo + Ao = 1. 





models 


ratio r 




OCDM (0.3,0.0) 


2.67 




LCDM (0.3,0.7) 
Inhom. « 2 = 2.2) 
Inhom. (C 2 = 2.0) 
Inhom. (C 2 = 1.8) 


25.0 


0.5 


141 




269 




586 




OCDM (0.3,0.0) 


19.5 




LCDM (0.3,0.7) 
Inhom. « 2 = 2.2) 
Inhom. (C 2 = 2.0) 
Inhom. (C 2 = 1.8) 


442 


0.4 


4.91 x 10 3 




1.33 x 10 4 




4.56 x 10 4 



3.2. Inhomogeneous cosmological models 

Here, we consider the inhomogeneous model with inner and outer homogeneous 
regions, as described in §1. For use of the Press-Schechter formalism in the two 
regions, we assume different spherical collapses in the two regions separately. Ac- 
cordingly, when we use Eq. ()2-15|) for S v , the parameters I?J and in the ho- 
mogeneous models are used in the inner and outer regions, respectively. Because 
the inhomogeneity is assumed to form from an adiabatic perturbation, the ratio of 
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.5 1 1.5 



Fig. 1. dN I (dzdO) for ct$(0)(2q = 0.5 in homogeneous models. L and O denote the models (LCDM 
and OCDM) with (J? , A ) = (0.3, 0.7) and (0.3, 0), respectively, in which r = 0.25 and h = 0.7. 
S and T denote the Einstein-de Sitter models, in which r = 0.5 and 0.25, respectively, and 
h = 0.5. 

the baryon density parameter to the total density parameter is equal in the two 
regions; that is, we have Q^/fi^ = Q^/Qq. This ratio has the value 0.116 for 
Q^hi 2 = 0.017. In this case, we have J?g/in 2 = 0.029, which is consistent with 
the value ~ 0.025 derived from the measurement of the deuterium abundance us- 
ing high-redshift QSOs, 40 ) while the inner value J?g/iT 2 = 0.017 is also consistent 
with the locally observed value. 41 ) Moreover, the value of p l vir corresponds to J7q 
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1 1 1 1 1 1 1 


1 


1 \ \ 




1 \ L \ 




-Is \ \ 

l^iis^ 1 1 1^ 1 


1 1 1 1 1 1 



1.5 



Fig. 2. dN / (dzdO) for ct$(0)(2q = 0.4 in homogeneous models. L and O denote the models (LCDM 
and OCDM) with (J? , A ) = (0.3, 0.7) and (0.3, 0), respectively, in which r = 0.25 and h = 0.7. 
S and T denote the Einstein-de Sitter models, in which _T = 0.5 and 0.25, respectively, and 
h = 0.5. 



for each i, and for the angular-diameter distance d\, we use the expression in the 
inhomogeneous model given in previous papers. 27 )' 29 ) 

The inner value of ag, [cs(O)] 1 , can be determined through local observations, in 
the same way as in homogeneous models using 



[<r 8 (0)] I (fto) P = °- 5 or °- 4 - 



(3-1) 
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.5 1 1.5 



Fig. 3. dN/(dzdO) for <7s(0)J?o = 0.5 in an inhomogeneous model (I) with £ = 2.0 in comparison 
with the two homogeneous models L and O. 

In the outer region II, on the other hand, it is to be noted that [o"8(0)] n is larger than 
the local value ([fs(0)]ioc) given in the inner region by the local observations, because 
the growth rate of density perturbations and the two-point correlation functions 
(of clusters) in the outer region are larger than those in the inner region, due to 
the difference between the model parameters. As a result, the present correlation 
function in the outer region also is larger than that in the inner region for a given same 
initial amplitude of density perturbations. Here, from the definitions of <Js(0) and 
the two-point correlation function £(R) (see Suto's review paper, 42 ) for instance), 
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Fig. 4. dN/(dzdfi) for <7g(0)f2Q = 0.4 in an inhomogeneous model (I) with £ = 2.0 in comparison 
with the two homogeneous models L and O. 



U8(0) and [C(-R)] 1 ^ 2 for a fixed distance R are proportional for a given same functional 
form of the power spectrum. For the model parameters we use in the regions I and 
II, the ratio of the correlation functions (= Q 2 ) is found to be of order 2, according 
to the theoretical analyses of Ref 32). For the local value of 08(0) specified by 

ks(0)]ioc(^ I ) p = 0.5 or 0.4, (3-2) 

therefore, we can adopt the value of in the region II given by 

M0)] H = CM0)]ioc, (3-3) 
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.5 1 1.5 2 



Fig. 5. dN / (dzdfi) for <r$(Q)Q% = 0.5 in inhomogeneous models with C/ 2 = 2.2,2.0 and 1.8, which 
are denoted by a, b and c, respectively. 

with C 2 ~ 2. Then we have M0)]7M0)] n ~ 0.3~ 0A5 /V2 = 1.22 for ( 2 = 2. 

In Figs. 3-6, the types of behavior of dN/ {dzdfi) for inhomogeneous models are 
compared with those for homogeneous models. Here we use units of deg -2 . In Figs. 3 
and 4 (for crg(0)J?Q = 0.5 and 0.4) we compare the case £ 2 = 2 with the homogeneous 
models LCDM and OCDM. It is found from Figs. 1-4 that for all z, the number 
density of clusters in the inhomogeneous model is much larger than that in SCDM 
and that for z < 0.25 (0.21), it is larger than that in LCDM for <7 8 (0)flg = 0.5 (0.4). 
Also, for z = 0.6, it is 1/2 (1/3) of that in LCDM for <7 8 (0)flg = 0.5 (0.4). Hence, 
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.5 1 1.5 



z 

Fig. 6. dN / (dzdfi) for u g (0)i?o = 0.4 in inhomogeneous models with ( 2 = 2.2,2.0 and 1.8, which 
are denoted by a, b and c, respectively. 

it is concluded that SCDM can be ruled out, because of the observed existence of 
clusters with z = 0.6 - 0.8. However, the present inhomogeneous model may be 
consistent with the observed existence of high-redshift clusters, like LCDM. 

In this inhomogeneous model, it is found on the other hand that for z ~ 0.1, the 
cluster number density in the outer region is larger by a factor of 1.5 - 2.0 than that 
in the inner region, and it is large also in comparison with the number density in 
LCDM. In Figs. 5 and 6 (for cr 8 (0)J?g = 0.5 and 0.4) the cases with ( = 2.2, 2.0 and 
1.8 are compared with the case of LCDM. It is found that the cluster number density 
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in the outer region depends strongly on the value of (, and that it decreases as a 
function of z more slowly for larger <ts(0). If C were equal to 1.0, the cluster number 
density in the outer region would reduce to that in the S model. The observation of 
the cluster number density for z ~ 0.1, therefore, may provide a stringent constraint 
on the value of £• 

The physical explanation of the behavior discussed above is basically that more 
clusters form in the outer region because the growth rate of density fluctuations 
there is larger than in the inner region and is characterized by a value of ( larger 
than unity. 

In Table I, values of the ratio r are listed for <78(0)I?q = 0.5 and 0.4 and for 
Q 2 = 2.2, 2.0 and 1.8. This shows that the ratio r increases with Q and that r in the 
inhomogeneous model with C/ 2 = 2.0 is larger by a factor ~ 10 than r in the LCDM. 

§4. Concluding remarks 

In the inhomogeneous model, the cluster number density in the outer region is 
larger than that in the inner region, as shown in the previous section. On the other 
hand, the galactic number densities in the two regions are nearly equal, as found in 
observational galactic surveys and studied theoretically in inhomogeneous models. 
We find that the number of galaxies within clusters in the outer region is smaller 
(by a factor of 1.5 2.0) than that in the inner region, and therefore the mass- 
luminosity ratio M/L in the outer region may be larger than that in the inner region 
by the same factor, as long as L of clusters is produced by galaxies within clusters. 
This result is interesting in connection with the observation of M/L of clusters made 
by Bahcall and Comerford 43 ) (cf. their Table I), in which we find that the mean 
value of M/L for clusters with z > 0.07 is larger by a factor of ~ 1.7 than that for 
clusters with z < 0.07. 

For the observed values of er 8 (0) , the values obtained from the cluster abundance 
in the region z < o.l 3 )' 35 )' 36 ) seem to be larger by a factor ~ 1.5 than the values in 
the region z > O.l. 38 )' 39 ) The difference between as(0) in these two regions may be 
due to the uncertainties, but if it is real, it may suggest the inhomogeneity of o"g(0). 
This is also interesting from the viewpoint of our inhomogeneous models, because 
o-g(O) in the inner region is larger by a factor ~ 1.2 than crg(0) in the outer region. 

The above-described characteristic behavior of the cluster abundance in the 
outer region of our inhomogeneous models will, moreover, be similarly found not 
only through the Sunyaev-Zeldovich effect, but also through X-ray and usual optical 
observations. A comparison of the theoretical and observational results for the latter 
methods with 2 > £ > 1 may give some interesting and severe constraints on the 
value of £• This must be studied subsequently. 

If the observed inhomogeneity of the Hubble constant is indeed real, it cannot 
be explained by homogeneous models but only by inhomogeneous models, while 
the accelerating behavior of SNla can be explained also by models with a Hubble 
constant inhomogeneity. Therefore, we finally add that more accurate measurements 
of the Hubble constant in gravitational lensing and SZE are desired. 44 ) The data for 
z = 1.7 26 ) are naturally consistent with the inhomogeneous model with (I? 1 ,, = 
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(0.3, 1.0), 30 - 1 while the consistency with LCDM has been argued as being due to the 
possible contribution of gravitational lensing to these data. 45 ^ 47 ^ To clarify this 
difference, we must wait for the day when many data of SNla with z > 1.5 are 
obtained. 
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